

enKnmis 


Professor J. V. Breakwell 
’ A. E. Bryson, Jr. (Coord 
Professor G. F. Franklin 


December 1971 


under 

Research Grant NASA-NgL-05-020-007 
from the 

National Aeronautics and Space Administration 

((NASA-CR-127268) GUIDANCE AND CONTROL OF N72-27679 

FLIGHT VEHICLES Annual Progress Report 
Jo Vo Breakwell, et al (Stanford UniVo) 

Deco 1971 18 p CSCL 17G Unclas 

! G3/21 15728 







ANNUAL PROGRESS REPORT 


December 1971 


on 


GUIDANCE AND CONTROL OF FLIGHT VEHICLES 


in the Departments of 

Aeronautics and Astronautics, and Electrical Engineering 


Principal Investigators 

Professor J. V. Breakwell 
Professor A. E. Bryson, Jr. (Coordinating) 
Professor G. F. Franklin 


under 


Research Grant NASA-NgL-05-020-007 


from the 


National Aeronautics and Space Administration 


This report summarizes progress during the past twelve months under 
continuing research grant for the period beginning January 1971. 



TABLE OF CONTENTS 


Page 


A. Studies in Attitude Control and Orbital Rendezvous 
(Supervised by J. V. Breakwell) ... 


1. Attitude Control of a Halo Orbit (J. V. Carroll) 

2. Minimum-Fuel Orbital Rendezvous between Neighboring 

Low Eccentricity Orbits (J. Jones) . . 

B. Studies in Optimization, Estimation, Identification, and 

Control (Supervised by A. E. Bryson) 

1. Estimation of the Local Attitude of Orbiting Spacecraft 

(W. Kortiim) 

2. Optimum Three Dimensional Turns for a Supersonic Air- 
craft (J. K. Hedrick) 

3. Square Root Filtering and Smoothing for Discrete Pro- 
cesses (P. G. Kaminski) 

» 

4. Combining VOR/DME Measurements with Air Data and with 

Inertial Data (J. C. Bobick) 

5. Optimum Inputs for Identification of Linear System 

Parameters (D. B. Reid) 

6. Guidance Logic for Orbit Transfer Using Low Thrust 

(L. J. Wood) 

7. Automatic Landing Logic for Aircraft (W. E. Holley) . . 

C. Studies in Identification (Supervised by G. F. Franklin) . . 

1. Identification of Linear System Parameters from EEG 

Signals (D. H. Utter) 

2. Identification and Modelling of Discrete, Stochastic 

Linear Systems (D. S. Spain) 

3. Equi control lability and the Model-Following Problem 

(R. T. Curran) 


1 

1 

3 

3 

3 

3 

4 
6 

6 

7 

8 

8 

14 

14 



RESEARCH RESULTS 


A. Studies in Attitude Control and Orbital Rendezvous 
Supervised by Professor J. V. Breakwell 

1. Attitude Control of a Halo Satellite (J. V. Carroll) 

This research, completed in 1971, examined the stability of a 
dual-spin halo satellite in the presence of radiation pressure torques 
as well as gravity-gradient torques. An almost spherical exterior 
surface is assumed, either rotating or non-rotating, while the inertia 
tensor is assumed to be axially symmetric. 

For a general non-symmetric non-spinning satellite, three-axis 
orientation may be controlled in the presence of radiation pressure 
torques by moving interior masses and hence also the location of the 
overall center of mass relative to the center of the exterior surface. 
Since no solar torque is obtained (to first order in the non-sphericity 
of the surface) about the sun direction, the time constant associated 
with such a control scheme is necessarily long. 

The final draft of doctoral dissertation is nearing completion. 

2. Minimum-Fuel Orbital Rendezvous between Neighboring Low Eccentri- 
city Orbits (J. Jones) 

- ✓ 

A computer program has been completed which can obtain, generally 
by iteration, the minimum-fuel N-impulse rendezvous (within a given 
duration). The equations are linearized about the target orbit, if 
circular, or about a neighboring circular orbit, otherwise. This 
linearization should provide very nearly optimal policies in many 
realistic situations. At the same time, the linearized problem has the 
great merit that all locally optimal rendezvous strategies, whether 
unique or not, are globally optimal. 

The iteration is started from the solution of a simpler rendezvous 
problem for which the solution is known. The simpler problems used 
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are of two types, (i) the "transfer" problem, in which the "phase in 
the final orbit is not specified, and (ii) the coplanar circle-to- 
circle rendezvous problem. In (i) the linearized analysis leads to 
three distinct types of strategy, one of which permits substantial 
freedom in the location of the impulses, and hence in the phase 
obtainable for the same fuel expenditure as in the transfer problem. 

If the desired phase is not within the range obtainable "for free,” 
the program enters the iterative phase, where careful step-size con- 
trol on the desired phase is combined with the possibility of increas- 
ing or decreasing the number of impulses in optimal fashion. For 
large desired phase change, starting type (ii) is preferable. As many 
as six impulses can be required, although five- and six-impulse situa- 
tions are sufficiently rare that four-impulse (or less) sub-optimal 

/ 

strategies may be computed with negligible fuel cost. 
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RESEARCH RESULTS 

B. Studies in Optimization, Estimation, Identification, and Control 

Supervised by Professor A. E. Bryson, Jr. 

1. Estimation of the Local Attitude of Orbiting Spacecraft (W. Kortum) 

This work, completed in 1970, has been published in Automatica , 
Vol. 7, pp. 163-180, 1971, under the title above. 

2. Optimum Three Dimensional Turns for a Supersonic Aircraft (J. K. 

Hedrick) 

This work, also completed in 1970, has resulted in three published 
papers: 

(i) "Minimum Time Turns for a Supersonic Airplane at Constant 
Attitude," Jour . Aircraft , Vol. 8, No. 3, March 1971, pp. 
182-187; 

(ii) "Three-Dimensional Minimum-Time Turns for a Supersonic 

Aircraft," Jour . Aircraft , Vol. 9, No. 2, February 1972, 
pp. 115-121; 

(iii) "Three Dimensional Minimum-Fuel Turns for a Supersonic 
Aircraft," Jour, i Aircraft , Vol. 9, No. 3, March 1972, 
pp. 223-229. 

3. Square Root Filtering and Smoothing for Discrete Processes (P. G. 

Kaminski) 

This work was completed in 1971 and published in SUDAAR Report 
No. 427, July 1971. A summary paper was published in IEEE Trans , on 
Auto. Control, Vol. AC-16, No. 6, December 1971, pp. 727-736, entitled 
"Discrete Square Root Filtering: A Survey of Current Techniques," by 

P. G. Kaminski, A. E. Bryson, Jr. , and S. F. Schmidt. 


The excellent numerical characteristics and reasonable computation 


4 


requirements of the square root approach make it a viable alternative 
to the conventional Kalman filter in many applications, particularly 
when computer word length is limited, or the estimation problem is 
badly conditioned. 

4. Combining VOR/DME Measurements with Air Data and with Inertial 
Data (J. C. Bobick — also supported by NSF traineeship and NASA 
Grant Ngr-05 -020-431) 

The first part of this work, dealing with air data, was presented 
at the AIAA Guidance, Control, and Flight Dynamics Conference, Hofstra 
University, August 1971, and is to be published in the Journal of Air— 
craft in June or July 1972 under the title "improved Navigation by 
Combining VOR/DME Information and Air Data. " A summary was given in 
last year's Progress Report. 

The second part of Mr. Bobick’ s research deals with updating or 
aligning an inertial navigation system in flight using VOR/DME signals, 
and is summarized below: 

The use of VOR/DME information with an INS has thus far been 
limited to the use of a DME/DME position fix to reset the position 
display (no filtering involved). This, of course, does not result in 
decreased velocity or platform attitude errors. 

Assuming a locally-level , north- pointing stable platform, we have 
studied the possibility of using the information from one or two VOR/DME 
stations as inputs to a filter which estimates the errors in the iner- 
tial navigation system. The states are the easterly and northerly com- 
ponents of the position and velocity errors, the angles of rotation 
about the east, north, and vertical axes which relate the desired and 
actual orientations of the platform, a drift in each of the three 
gyros, an error in the east and north accelerometers, and a bias error 
in each VOR and each DME. 

The estimates of the states determined by the filter can be used 
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to torque the platform to its desired orientation, or, as long as the 
platform attitude errors remain small, we can simply compensate for 
them without realigning the platform. Using such a filter, we have 
determined the RMS errors in the estimates of position, velocity, and 
platform attitude (cp^ = tilt about east axis, cp N = tilt about north 
axis, and cp^ = azimuth error) for a 30-minute flight using (a) one 
VOR/DME to update the INS, and (b) two DME's, to update the INS. The 
present method of ground alignment takes about 15 minutes with the 
final platform tilts, (cp^,)^, (cp^f > equal to about 0.005 degrees 

and the final azimuth error, (cp^)^ * equal to about 0.05 degrees. We 
conclude that the accuracy of in-flight alignment using two DME's is 
about the same as for ground alignment, whereas the accuracy of in- 
flight alignment using one VOR/DME is worse than ground alignment by a 
factor of 2 or 3. 

Although the need for initial in-flight alignment of inertial 
systems onboard commercial aircraft is questionable, the in-flight re- 
alignment of the system before leaving the U. S. airspace on a trans- 
oceanic flight (e.g. , from Los Angeles to London or Chicago to Hawaii) 
could result in significant improvements in navigational accuracy. 

Also, realignment of the system upon entering the U. S. airspace after 
a transoceanic flight could prove useful since the result would be ac- 
curate position, velocity, and attitude information in the terminal 
area. 

We found that the suboptimal filter resulting when the gyro 
drifts, accelerometer errors, and DME biases are neglected, performed 
nearly optimally during the periods of time the filter operates (20 to 
30 minutes). Hence, the results presented can be obtained by using an 
eighth order filter when using one VOR/DME or a seventh order filter 
when using two DME's to update an INS. 
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5. Optimum Inputs for Identification of Linear System Parameters 
(D. B. Reid — also supported by NASA Grant NgL-05-020-526) 

A major portion of this study was completed in 1971. This included 
determining a suitable criterion for evaluating inputs and numerical 
solution of two aircraft identification problems. The criterion used 
was to maximize the trace of the information matrix for the parameters 
to be identified. The problems solved were (a) optimum elevator input 
to identify the pitching moment damping derivatives and 

from pitch rate measurements, and (b) optimum aileron/rudder inputs to 
identify the yawing and rolling moment damping derivatives C nr> C np ,Cg r , 
Cg from yaw rate and roll rate measurements. All of these inputs 
were of the bang-bang type, so solutions required finding the switch- 
ing times from maximum to minimum deflection, or vice-versa. A report 
is in preparation. 

6. Guidance Logic for Orbit Transfer Using Low Thrust (L. J. Wood) 

This study was nearly completed in 1971. The major new result is 
a guidance scheme for non-autonomous systems with open end time, with 
a detailed example. 

Minimum-time low thrust orbit transfers were computed for a typical 
spacecraft from Earth to Mars and from Earth to Jupiter. Neighboring 
optimum feedback gains were determined for the minimum-time Mars trans- 
fer including a time-to-go logic for entering the gain tables and the 
nominal trajectory tables. This latter scheme was simulated on the 
computer and worked well even for very large initial errors (i.e. , 
deviations from the optimum nominal path). To our knowledge, this is 
the first demonstration of a workable minimum-time guidance scheme for 
a non-linear system of this complexity. 

An explanation was found for the unexpected retro-thrust period 
that occurs in the middle of interplanetary transfer path: low thrust 

achieves outward radial velocity partially by a radial component of 



thrust, but also by a tangential component to increase speed so that 
centrifugal force exceeds the sun's gravitational force. By the same 
token, to decrease radial velocity, a backward component of tangential 
thrust (retro-thrust) is used in the middle of the transfer path. 

A report is in preparation. 

7. Automatic Landing Logic for Aircraft (W. E. Holley) 

This investigation has only begun. Preliminary work indicates 
that the use of integral control concepts and/or filters to estimate 
the wind velocities combined with quadratic synthesis techniques pro- 
vides significant improvements in autoland accuracy in the presence of 
winds. Statistical wind models were studied in relation to the data 


available. 
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RESEARCH RESULTS 

C. Studies in Identification 
Supervised by Professor 
G. F. Franklin 

1. i dentification of Linear System Parameters from EEG Signals (D. H. 

Utter) 

The problem under investigation is to see if it is possible to 
identify parameters of a linear model of low order which will adequately 
describe the changes in the electroencephalograph process during anes- 
thesia. The results of the investigation can be used in on-line identi- 
fication of a parametric model for any stochastic process when only 
the output signal can be measured. We want to know what the important 
parameters are in terms of sampling rates, estimate accuracy, and the 
assumption of stationarity in making a low order model of a very complex 
process. / 

The particular problem arises from a study of the level of anes- 
thesia in patients undergoing surgery. It has been observed that 
electroencephalograph (EEG) recordings of these patients usually 
exhibit the characteristics of a stationary stochastic process over 
short time periods, and their spectra show recognizable changes with 
changes in anesthesia conditions. One standard procedure for EEG 
monitoring has been to process the signals by spectral analysis. This 
method preserves all the information in the signal but produces a 
large number of parameters: the signal power at many discrete fre- 

quencies. So much information would be overwhelming to the busy anes- 
thesiologist during surgical operations and very difficult to use in a 
control algorithm, should such be desirable. For these applications, 
a more meaningful procedure appears to be to represent the stochastic 
process by a parametric model consisting of a linear filter of low 
order driven by white noise with the EEG signal as the output. 

For the identification, it is assumed that the white noise is 
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normally distributed with unity variance. The filter output signal is 
measured at discrete time intervals of T seconds separation and 
stored as a finite-length sample of N data points. The parametric 
model is a linear difference equation representing the output at the 
i th time point as given by Eq. (1): 

n 

i Z- J 1- J 

j=l 

where u is the white noise input at the i th time point. The coef- 
x 

ficients a and b. are the parameters to be identified. The index 
j J 

n is the assumed order of both the denominator and the numerator of 
the discrete filter function and is a parameter to be chosen in the 
identification algorithm. 


b .u . 
J i~J 


j=0 


( 1 ) 


Results to date indicate that the a . coefficients can be es- 

J 

timated efficiently and quickly from an algorithm similar to one dis- 
cussed by Zetterberg [1]. The estimate is asymptotically unbiased. 

It is derived using the fact that the expected value of the product of 
future noise input and present filter output is zero. If we multiply 
(1) by y Q and take expectations, then for i > n the terms of u fc 
vanish. Thus the equations for i = n+l,...,m can be arranged in 
matrix form as Eq. (2): 


Here 


r = -Ra . 

is the n-vector of estimated 


( 2 ) 


a. coefficients, r 
J 


is a vec- 


tor of (m-n) autocorrelation coefficients such that r = (r , ,r „, 

— n+l n+z 
T 

...,r ) with m a parameter > 2n to be chosen later. The matrix 
m — 

R is comprised of (m-n) rows and (n) columns of autocorrelation 
coefficients such that each row of Eq. (2) forms the Eq. (3): 


r 


k 


r. . a. 
k-i i 


i=l 


(3) 
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Equations (3) are the well known Yule-Walker equations for the given 
assumptions [2], The autocorrelation coefficients, r are estimated 

K 

from the sample of N points by Eq. (4): 

N-k 

I (y i-y )<y i +k -y ) (4) 

i=l 

where y is the sample mean. The constant 1/N is used rather than 
1/N-k , the constant used in the standard formula for computing auto- 
correlation coefficients from a finite-length sample, and is equivalent 
to using a weighting factor of N-k/N in order to give less importance 
to the correlations of greatest time shifts. The formula may also be 
seen as a "window" on the data for estimation of r k . Using (4) in 
(2), the equations are solved by standard least squares methods to 
give the characteristic equation coefficients of the model. 

The a are converted into z-plane (sampled data) roots z ± by 
solving for the complex roots of the polynomial equation 

n 

1 + a^z i = 0 . (5) 

i=l 


Finally, in order to relate the parameters to known spectral features, 
the equivalent s-plane (continuous) roots of the denominator are com- 
puted by Eq. (6) for real roots and Eqs. (7) and (8) for complex 
roots: 


s. = - in z, , z . real 

l T i i 


Re(s i ) = J £n l z l 

Im(z . ) 
i 

Re(z. ) 


( 6 ) 

(7) 


Im(s i ) = i arctan 


( 8 ) 
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This algorithm for estimating the a^ coefficients has been applied 
to steady-state EEG data from different patients at various levels of 
anesthesia, as identified by Dr. J. Bimar, an anesthesiologist who 
specializes in EEG monitoring. The data samples contained N = 1024 
data points at intervals of T = 1/64 seconds. Several values of the 
parameters n and m were tried. Visually consistent results were 
obtained for n = 5 and m = 200 . Figure 1 shows the s-plane root 
locations obtained from 17 samples from a single record of steady- 
state data at a moderate anesthesia level designated by Dr. Bimar as 
MAC1. Typically, three of the five roots identified are plotted for 
each sample: one real root and two upper-half -plane complex roots. 

Figure 2 illustrates the s-plane root locations identified from 16 
samples of another record at a deeper anesthesia level designated MAC2 
by Dr. Bimar. In each figure the lower of the two complex roots is 
plotted as "+" while the real root and the higher frequency complex 
root are both indicated by "x". The fact that these roots cluster is 
taken as a measure of the stability of the identification technique. 
The fact that the clusters are in different parts of the s-plane is 
taken as a sign that these parameters can be used to distinguish be- 
tween these different states of anesthesia. 

The investigation continues into effective methods for estimating 

the b. coefficients. 

J 
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2. Identification and Modelling of Discrete, Stochastic Linear Systems 
(D. S. Spain) 

This research was completed in 1971 and published with the above 
title as SEL Technical Report No. 6302-10. The basic problem dealt 
with is the identification of stochastic, discrete, multivariable 
linear systems from input -output data. The problem definition includes 
the possibility of both deterministic and stochastic inputs, although 
all the noise sources are restricted to be gaussian and white. Using 
an innovations formulation of the maximum likelihood criterion, proba- 
bility one convergence of the impulse response matrices to their true 
values can be obtained. From these matrices, a canonical form is de- 
veloped for multivariable linear systems which requires no structural 
information or other prior knowledge of the system (although the re- 
sults are derived in such a way that such knowledge can certainly be 
used if it is available). This canonical form is also useful in the 
least squares identification of multivariable systems. Two very satis- 
factory methods of identifying the dimension of a system are presented, 
one based on the whiteness of the resulting error process, and the other 
based on the relative decrease of the cost function as the dimension of 
the model increases. Conditions are also given which the deterministic 
input must satisfy to guarantee the probability one convergence of the 
impulse response matrices, and easy methods are pointed out of con- 
structing input sequences which satisfy these conditions. Finally, 
both a working program to perform the identification, and suggestions, 
based on computational experience, for possible improvements, are 
presented. 

3 . Equicontrollability and the Model -Following Problem (R. T. Curran) 

This research was completed in 1971 and published as SEL Technical 
Report No. 6303-2 under the title listed above. It introduces the idea 
of equicontrollability and studies its application to the linear time- 
invariant model-following problem. The problem is presented in the 
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form of two systems: generically called the plant and the model. The 

requirement is to find a controller to apply to the plant so that the 
resultant compensated plant behaves, in an input -output sense, the 
same as the model. All systems are assumed to be linear and time- 
invariant. 

The basic approach used is to find suitable equicontrollable 
realizations of the plant and model and to utilize feedback so as to 
produce a controller of minimal state dimension. The concept of equi- 
cont reliability (introduced here) is a generalization of control 
canonical (phase variable) form applied to multivariable systems. It 
allows one to visualize clearly the effects of feedback and to pin- 
point the parameters of a multivariable system which are invariant 
under feedback. 

The basic contributions contained in this work are: (1) the 

development of equicontrollable form, (2) solution of the model-fol- 
lowing problem in an entirely algorithmic way, suitable for computer 
programming, and (3) resolution of some questions on system decoupling, 
along with the application of the above algorithm to accomplish decoupling. 



